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Abstract 

The inclusive double-spin asymmetry, A^, has been measured at COMPASS in deep- 
inelastic polarised muon scattering off a large polarised NH3 target. The data, collected 
in the year 2007, cover the range > 1 (GeV/c)^, 0.004 < x < 0.7 and improve the 
statistical precision of gi{x) by a factor of two in the region x < 0.02. The new pro- 
ton asymmetries are combined with those previously published for the deuteron to extract 
the non-singlet spin-dependent structure function g^^{x,Q'^). The isovector quark den- 
sity, Aq^lx, Q"^), is evaluated from a NLO QCD fit of . The first moment of Aga is 
in good agreement with the value predicted by the Bjorken sum rule and corresponds to a 
ratio of the axial and vector coupling constants \gA/gv\ = 1.28±0.07(stat.)ib0.10(syst.). 
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In previous publications the COMPASS collaboration has presented new accurate values 
of the longitudinal spin asymmetry of the deuteron, Af, covering a large range of x (0.004 < 
X < 0.7) in the region of deep inelastic scattering (DIS) > 1 (GeV/c)^ |[Dl2l- These new 
values have led to an improved determination of the spin structure function gf{x, Q"^) in the 
low X region where only the SMC measurements existed before [3J. The first moment of gf{x) 
has also provided a more accurate value for the matrix element of the flavour singlet axial 
current = 0.35 ± 0.03(stat) ± 0.05(syst.) at = 3 (GeV/c)^, confirming the rather small 
contribution of the quark spins to the nucleon spin. 

In this letter we present new COMPASS results for the inclusive double-spin asymmetry 
of the proton, and for the spin structure function g\{x, Q"^) measured in the same kinematic 
range. In combination with the deuteron data, these results yield an evaluation of the isovector 
quark density, /^qz{x, Q"^) = Au — Ad, and its first moment which in turn provides a test of the 
Bjorken sum rule. 

The proton data were collected in 2007 using a three target cell configuration and the 
upgraded COMPASS spectrometer, as described in [4J. The polarised target material is the am- 
monia previously used in the SMC experiment Q. The polarisation was about 90% in absolute 
value, measured with a relative error of ± 2% [6J. The muon beam has a natural polarisa- 
tion of about —80%. The energy of the incoming muons is constrained to be in the interval 
140 < < 180 GeV and their polarisation is known with a relative precision of ± 5%. All 
events used in the present analysis are required to have a reconstructed primary interaction ver- 
tex (defined by the incoming and the scattered muon) inside one of the target cells. In order to 
cancel out the muon flux normalisation in the asymmetry calculation, incident muons are only 
accepted when their extrapolated trajectory crosses all three cells. 

For most events the trigger is based on a combination of hodoscope signals defining the 
trajectory of the scattered muon. In addition to these "inclusive triggers", low x events are also 
selected by an additional condition on the energy deposit in the hadron calorimeter, which is 
then used as a "semi-inclusive trigger". At large x and Q"^ most events are selected by conditions 
on the calorimeter signal only, without any input from hodoscopes. For this "calorimeter-only 
trigger" as well as for the semi-inclusive one, the presence of a reconstructed hadron trajectory 
compatible with the calorimeter information is required. 

The kinematic region is defined by requiring the photon virtuality > 1 (GeV/c)^ and 
the fractional energy y transfered from the beam muon to the virtual photon to be between 0.1 
and 0.9. The region which is most affected by radiative corrections is eliminated by the cut 
y < 0.9. The total sample after all cuts amounts to 85.3 million events. 

The longitudinal virtual-photon proton asymmetry, A\, is evaluated from the numbers 
of events collected in the different target cells by the method used in our previous analyses 
of deuteron data [2113 • Neighbouring target cells are polarised in opposite directions and data 
from both target spin orientations are thus recorded simultaneously. The lengths of the cells 
are chosen so that the two samples collected with opposite spin orientations have in average 
the same acceptance, which limits the risk of false asymmetries. The target spin directions are 
reversed once per day by rotating the magnetic field and a few times per year by changing the 
microwave frequencies used for dynamic nuclear polarisation. The asymmetries are calculated 
from the numbers of events in cells with opposite spin orientations collected before and after a 
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field rotation so that flux and acceptance factors cancel out. 

Radiative corrections are applied separately to the asymmetries obtained for the inclu- 
sive triggers and to those obtained for the semi-inclusive and calorimeter-only triggers because 
radiative elastic events contribute only to the former ones. Another correction is applied to 
account for the polarisation of the ^''N nucleus. For this spin 1 object the correction is propor- 
tional to Af{x) and affected by factors accounting for the number of ^^N nuclei vs. H atoms, 
the alignment of the proton spin vs. the ^^N spin and the ratio of ^"^N to H polarisations [|71: 
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The corrections for the various intervals of x are given in the appendix. They are of the order 
of 0.01 for X > 0.35 and are mainly important for the evaluation of the first moment of the spin 
structure function rf(Q^) = g{{x,Q'^)dx. 

The target dilution factor is given by the ratio of the cross-section for the polarisable 
protons to that of all nuclei in a target cell. The values for the NH3 target are shown in Fig. [T] as 
a function of x, for inclusive and hadron triggers. They are about 14% in the medium x region, 
with a rise at large x due to the reduced cross section on heavy targets in this region, and a drop 
at low X for inclusive triggers due to the contribution of radiative elastic events on the proton. 
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Figure 1 : The dilution factor / of the NH3 polarised target as a function of x for inclusive and hadron 
triggers [3]. The values of f{x) are averaged over the range of the corresponding triggers. 

The new values of A^{x) are shown in Fig. [2] in comparison with results from previous 
experiments. The of different points at any fixed value of x varies considerably since the in- 
cident energy of the various experiments ranges from 6 to 200 GeV. The fact that all results align 
reasonably well on a single curve illustrates the well known observation that the dependence 
of is very weak in the DIS region. This is further illustrated in Fig.[3]which shows A^(x, Q"^) 
as a function of for the COMPASS data. No significant dependence is observed in any 
interval of x. 

The systematic errors of the COMPASS results for A^ are shown by the band at the 
bottom of Fig. [2] and listed in Table 1. They contain the contributions due to the uncertainties 
on the target polarisation, the beam polarisation, the dilution factor and the ratio R = ai/o-T 
jH used in the depolarisation factor, which are equal to 2, 5, 1 and at most 3%, respectively. 
Combined in quadrature, these uncertainties amount to a systematic error of at most 6% of 
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Figure 2: The asymmetry A^{x) as measured by COMPASS and previous results from EMC [9|, SMC 
Dl, HERMES [10], SLAC E143 [11], E155 [12] and CLAS [13] at > i (GeV/c)^. The cut W > 
2.5 GeV has been applied to select DIS events in the CLAS data. Only statistical errors are shown with the 
data points. The band at the bottom shows the estimated size of the systematic errors for the COMPASS 
data. 



the quoted value. The error due to the neglect of the transverse asymmetry, Af, is less than 
0.002 in the full range of x. Another possible contribution to the systematic error is due to 
false asymmetries generated by instabilities in some components of the spectrometer. Such 
effects have been searched for in faked configurations, where the physics asymmetry does not 
contribute, and were found to be compatible with zero. The comparison of results obtained 
with opposite orientations of the target field also does not show any significant difference. The 
possible error due to false asymmetries has been estimated by a statistical test performed on 
the distribution of asymmetries extracted from 46 subsamples. Time-dependent effects which 
would lead to a broadening of these distributions were not observed. As a consequence the limit 
Csyst < 0-'i7o'stat was obtained at the level of one standard deviation. The different contributions 
to the systematic error are summarised in Table [2l 

The longitudinal spin structure function of the proton is obtained from by the 
relation 

where F| is the spin independent structure function. The values obtained with the SMC pa- 
rameterisation of the world data on [3] and the parameterisation of R already used in the 
depolarisation factor are listed in Table [H with their statistical and systematic errors. They are 
also shown in Fig.|4]in comparison with the SMC values. It can be seen that the COMPASS data 
improve the statistical precision at least by a factor of two in the low x region, covered only by 
the two experiments shown here. The new data points are compatible with a constant gi{x) for 
0.004 < X < 0.04 and do not show evidence either for an increase or a decrease when x — )■ 0. 
This observation remains valid when the data points are moved to a common Q"^ according to the 
fits quoted in Ref. [|2]| and the constant value is found to be = 0.48±0.03(stat.) ±0.04(syst.) 
atg2 = 3(GeV/c)2. 
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Figure 3: Values of as a function of in intervals of x. The errors are statistical only. The solid 
lines show the results of fits to a constant. 



The non-singlet spin structure function 

g'^\xM'')=9\{x.Q^)-g1{x.Q^) (3) 

is of special interest because its dependence is decoupled from the singlet and the gluon spin 
densities: 

g^M') = i£^C^^(^,a.(g^))Ag3(a:',Q^) (4) 

where C^^ is a Wilson coefficient function and Ags the isovector spin density. Consequently a 
fit of the Q'^ evolution of gf^ requires only a small number of parameters to describe the shape 
of Ag3(x) at some reference . According to the Bjorken sum rule the integral of g^^ at any 
fixed is proportional to the ratio gA/gv of the axial and vector coupling constants and given 
by the relation 



rf^(g2) = 1 i££L c^'{Q'') (5) 
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where the non-singlet coefficient function Cf^'^(Q^) has been calculated in perturbative QCD 
up to the third order in asiQ"^) [14]. The comparison of the value of \gA/9v\ obtained from the 
data with the one derived from neutron (3 decay (\gA/9v\ = 1.2694±0.0028 El) thus provides 
a test of the Bjorken sum rule, free of systematic errors arising from uncertainties on the gluon 
helicity distribution. 
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Figure 4: The spin structure function g^{x,Q'^) vs. x as measured by COMPASS at the of each 
measured point. Previous results from SMC lO are shown for comparison. The errors are statistical only. 



In the present analysis, the values of g^^ are obtained as 



1 - l.btu 



D 



(6) 



where the values of gf are taken from Ref. [2J and the deuteron D-state probability, lod, is 

0.05 ±0.01 dm. 

The proton and deuteron data have been obtained at the same {x, Q"^) points, which avoids 
the need of interpolation. A correction has however to be applied to the deuteron data for the 
admixtures of and in the ^LiD target material, which was not taken into account in Ref. 
[[2J. The ratios of isotopes ''Li/^Li and H/D were found to be 4.4 % and 0.5 % respectively IfTTl 
and ^Li and are both polarised at more than 90% [18J. The resulting corrections to are 
given in the appendix. They are negligible at low x but reach 0.015 at x = 0.50 and reduce the 
first moment Vf by 0.002. 

In the present analysis, = 3 (GeV/c)^ has been taken as reference Q'^ and the follow- 
ing parameterisation has been used for Aga: 

Ag3(x) = n^-i \ . (7) 

As for our previous analysis [2] the QCD fit at NLO of the COMRASS values of g^^ to Eq.(4) 
has been performed with two different programs, the first one working in the {x,Q'^) space [fT9l , 
the second one in the space of moments [|20|. Both programs give the same values of the fitted 
parameters and similar -probabilities. The fitted parameters obtained are listed in Table [3l 
The exponent 0:3 is in the range of Regge pole predictions [21] and the integral rj^ is in excellent 
agreement with the Bjorken sum rule prediction \gA/ gv\- The fitted distribution of xg^^{x) and 
the data points moved to the reference are shown in Fig. [5] (left). 

The integral of g^^ has also been evaluated from the measured values in the range 
0.004 < X < 0.7 with additional low and high x contributions taken from the fit (Table ID). 
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Table 1: Values of and g'^ as a function of x with the corresponding average value of Q^. The first 
error is statistical, the second one systematical. 



Beam polarisation 




0.04/0.8 = 5.0% 


Target polarisation 




2% 


Depolarisation factor 


AD{R)/D{R) 


2.0-3.0% 


Dilution factor 




1 % 


Total 




~ 0.06^1 


Transverse asymmetry 


7]-A2 


< 2.0 X 10-3 


Rad. corrections 


AAf^ 


10-^ - 10-3 


False asymmetry 


Afalse 


< 0.47 ■ AAf 



Table 2: Decomposition of the systematic error of yl^ into multiplicative (top) and additive (bottom) 
contributions. 

It is observed that about 92% of the first moment T^^ comes from the measured region. The 
dependence of the first moment of g^^ on its lower limit is shown in Fig. [5] (right). As already 
observed in the HERMES analysis [fTOl . the integral does not saturate at x ~ 0.01 — 0.02 while 
the value obtained at the lowest x accessible in the present analysis (0.180 ± 0.009) is less than 
one standard deviation below the value expected from the Bjorken sum rule (0.188). The value 
of \gA/ gv \ derived from the value of V^^ by Eq.(|5]) is identical to the one obtained from the fit 
and confirms the validity of the Bjorken sum rule with a statistical precision of 5%. 

The dominant systematic error on this result is due to the uncertainty of 5% on the beam 
polarisation, which is common to the proton and deuteron data and therefore translates directly 
into a 5% error on \gA/ gv\- Other contributions due to the target polarisation and the dilution 
factor are estimated to be ±0.04 and ±0.06 for the proton and deuteron terms, respectively. The 
resulting systematic error is ±0.10. The errors related to the fit or to the evolution of the data 
to a common are found to be negligible. In particular, it was checked that the same value of 
Qa/ gv is obtained when the reference is 1.0, 3.0 or 10.0 (GeV/c)^ although the exponent 
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Par am. 


Value 


^3 


1.28 ±0.07 


as 


-0.22 ±0.07 




o o_(_0.5 


xVnd 


14.4/12 


Prob. 


0.27 



X range 




- 0.004 
0.004 - 0.7 
0.7 - 1.0 


0.0098 

0.175 ±0.009 ±0.015 
0.0048 


- 1 


0.190 ±0.009 ±0.015 



Table 3: Results of the fits of Aq^^x) at 
g2 = 3 (GeV/c)2. 



Table 4: First moment rf ^ at = 3 (GeV/c)^ from the 
COMPASS data points. The contributions from the unmea- 
sured regions are estimated from the NLO fit to Qi^; their 
errors are negligible. 




g^S(x)dx 


™"0.2 
0.18 




0.16 




0.14 




0.12 




0.1 




0.08 




0.06 




0.04 




0.02 










-Bjorken sum rule 



COMPASS data 



10" 



10' 



X 1 

min 



Figure 5: Left: Values of gi^^{x) at Q"^ = 3(GeV/c)^, derived from the COMPASS measurements 
of and Af and result of a three parameter QCD fit at NLO. The errors are statistical only. Right: 
as a function of Xmin as obtained from the COMPASS data points. The open circle at 
X = 0.7 is obtained from the fit. The arrow on the left side shows the value expected for the full range 
< X < 1 with \gA/gv\ = 1-269 ifBH . 



as varies from —0.15 to —0.28 when Q'^ is moved from 1 to 10 and the shape of g^^{x) thus 
becomes quite different. 

The test of the Bjorken sum rule performed in the present analysis of the COMPASS 
proton and deuteron data is thus mainly limited by systematics: 

\gA/gv\ = 1-28 ± 0.07(stat.) ± 0.10(syst.), (8) 

to be compared with the value 1.2694 ± 0.0028 derived from neutron f3 decay [fTSl . 
The COMPASS value of T^^is in good agreement with the one obtained by the SMC (0.198 ± 
0.023 at = 10 (GeV/c)^) [[3]| and improves the statistical precision by a factor 2.5. The 
cumulative integral of g^^ truncated at Xmin = 0.021 is equal to 0.1583 ± 0.0085(stat.) ± 
0.014(syst.), in good agreement with the HERMES value of 0.1484±0.0055(stat.)±0.016(syst.) 
obtained at = 5 (GeV/c)^ ifTOl. 

In conclusion, the COMPASS collaboration has performed new measurements of the lon- 
gitudinal spin asymmetry of the proton, covering a large range of x (0.004 < a; < 0.7) in the 
DIS region, > l(GeV/c)^. The new data improve the statistical precision in the low x region 
by a factor of 2-3 and show no evidence either for an increase or a decrease of the spin structure 
function g^ in this region. In combination with the previously published results on the deuteron, 
the new data improve the evaluation of the non-singlet spin structure function g^^ and provide 
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a test of the Bjorken sum rule, which is satisfied within one standard deviation of the statistical 
uncertainty. 
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Appendix 



X range 


Corr. to 


Corr. to Af 


0.004 


- 0.005 


0.000 


0.000 


0.005 


- 0.006 


0.000 


0.000 


0.006 


- 0.008 


0.000 


0.001 


0.008 


- 0.010 


0.000 


0.001 


0.010 


- 0.020 


0.000 


0.001 


0.020 


- 0.030 


0.000 


0.001 


0.030 


- 0.040 


-0.000 


0.002 


0.040 


- 0.060 


-0.001 


0.002 


0.060 


-0.100 


-0.001 


0.003 


0.100 


-0.150 


-0.003 


0.004 


0.150 


- 0.200 


-0.004 


0.006 


0.200 


- 0.250 


-0.005 


0.008 


0.250 


- 0.350 


-0.006 


0.010 


0.350 


- 0.500 


-0.009 


0.013 


0.500 


- 0.700 


-0.014 


0.017 



Table 5: Corrections to the COMPASS spin asymmetries and Af due to the the ^^N polarisation 
and to the admixture of ^Li and into the ^LiD target material. In both cases the correction must 
be subtracted from the measured asymmetries. The corrections to A^ are already applied to the values 
quoted in the present letter. 
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